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EXECUTIVE  SUMMARY 


The  fourth-order  spectrum  represents  a  new  way  of  extracting  information  from  data.  Used  in  conjunc¬ 
tion  with  the  second-ora  >r  spectrum  (or  power  spectrum)  the  fourth-order  spectrum  reveals  otherwise  hidden 
relationships  that  are  important  for  classification  of  transient  signals  in  passive  sonar  and  for  classification  of 
resonances  in  active  sonar.  In  addition,  the  fourth-order  spectrum  can  extract  range  and  Doppler  information 
from  an  active  sonar  return  under  conditions  which  render  the  conventional  spectrum  useless.  Moreover,  the 
potential  for  signal-to-noise  ratio  improvements  in  Gaussian  noise  environments  is  high.  Theoretically,  the 
fourth-order  spectrum  eliminates  all  additive  Gaussian  noise  from  its  final  result.  The  practicality  of  this 
theoretical  fact  is  currently  being  considered. 

The  drawbacks  in  using  the  fourth-order  spectrum  in  practice  are:  1)  it  requires  more  processing 
capacity,  and  2)  it  will  require  training  in  order  to  properly  interpret  the  results,  unless  the  detection 
and  classification  functions  are  automated.  However,  these  drawbacks  have  been  minimised  in  this  report 
by  the  utilisation  of  a  special  case  of  the  fourth-order  spectrum.  Nevertheless,  all  the  advantages  of  the 
fourth-order  spectrum  discussed  above  have  been  maintained. 

The  following  points  are  discussed  in  more  detail  in  the  report: 

1)  A  new  class  of  multidimensional  non-Gaussian  density  functions  are  introduced.  This  class  represents 
physically  meaningful  signals,  since,  independent  data  are  not  required.  As  a  special  case,  the  class  reduces 
to  a  multidimensional  Gaussian  density,  so  a  wide  variety  of  density  functions  can  be  represented.  In  the 
report  fourth-order  spectra  were  obtained  for  the  class  of  multidimensional  non-Gaussian  density  functions. 
In  this  way,  it  was  shown  that  the  fourth-order  spectrum  could  differentiate  between  a  Gaussian  and  a  non- 
Gaussian  process.  This  is  important  because  many  noise  sources  in  the  ocean  are  Gaussian.  For  example, 
ambient  noise,  reverberation,  flow  noise,  ect.,  are  Gaussian.  Whereas,  signals  are  usually  non-Gaussian.  For 
example,  sinusoids,  transients,  active  sonar  transmissions,  ect.,  are  all  non-Gaussian. 

2)  An  active  sonar  classification  technique  relies  on  backscatter  resonances  to  differentiate  target  types. 
But  resonances  may  not  give  an  unique  classification,  since,  many  waveforms  could  have  the  same  resonances. 
In  order  to  show  the  potential  of  the  fourth-order  spectrum  for  classification,  an  idealised  example  of  creeping 
wave  backscatter  is  discussed.  In  this  example,  a  sum  of  two  sinusoids  and  a  mixture  of  two  sinusoids  have 
the  same  spectrum  (resonances).  However,  only  the  mixture  represents  the  physical  phenomenon  of  creeping 
wave  backscatter.  It  is  shown  that  the  fourth-order  spectrum  differentiates  between  a  sum  and  a  mixture  of 
sinusoids,  and  therefore,  can  aid  in  active  sonar  classification. 

3)  When  long  pulse  trains  are  transmitted  the  return  echo  is  usually  amplitude  modulated  from  target 
dynamics,  Doppler  spreading,  or  medium  effects.  Under  these  conditions,  the  spectrum  can  be  severly 
distorted  making  detection  and  classification  impossible.  Using  an  amplitude  modulated  coded  pulse  train, 
it  is  shown  that  the  fourth-order  spectrum  can  extract  range  and  Doppler  information  while  for  the  same 
conditions  the  spectrum  is  useless. 
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Fourth-Order  Spectra  of  Mixture  and  Modulated  Processes 


1.  Introduction 


Certain  underwater  acoustic  signals  can  be  modeled  as  mixture  processes.  For  example,  a  scattered  wave 
from  an  extended  target  will  often  consist  of  several  subwaves  with  different  characteristics.  Determining 
whether  or  not  a  received  waveform  is  from  an  extended  target  is  an  important  problem[l,2).  Another 
example  of  a  mixture  process  is  an  intermitteney  produced  by  a  complex  mechanical  source[3|.  Also, 
environmental  conditions  can  cause  several  signals  with  different  doppler  frequencies  to  be  mixed  together 
in  the  received  waveform.  In  all  of  these  examples  the  fundamental  properties  of  the  waveform  can  be 
characterised  by  a  mixture  process.  Other  signals  are  better  modeled  as  modulated  processes.  Wave 
propagating  in  random  medium  and  Doppler  spread  targets  are  two  examples. 

The  objectives  of  the  report  are:  1)  to  develop  a  parameterised  mixture  process  model,  similar  but 
more  general  than  reference  4,  that  characterises  the  phenomena  of  interest  and  extract  information  from 
the  model  by  estimating  Us  spectra  and  higher-  order  spectra  j5,6,7);  and  2)  to  estimate  spectra  and  higher- 
order  spectra  from  amplitude  modulated  processes. 

It  will  be  shown  in  the  report  that  the  mixture  process  is  in  general  a  non-Gaussian  process.  Therefore, 
higher-order  spectra  would  contain  additional  information  about  the  non-Gaussian  waveform.  In  our  initial 
application,  a  measurable  higher-order  spectrum  would  be  enough  to  classify  the  phenomenon  as  belonging 
to  a  mixture  process.  Specifically,  it  will  be  shown  that  a  special  case  of  the  trispectrum  [8,9]  can  determine 
whether  or  not  a  mixture  process  is  present,  whereas,  the  second-order  spectrum  cannot.  For  this  special 
case,  an  interesting  result  for  sinusoidal  mixture  processes  is  presented  and  demonstrated  by  simulation. 
For  the  modulated  case,  examples  are  presented  which  show  the  usefulness  of  the  fourth-order  cumulant 
spectrum  in  extracting  target  information  from  an  amplitude  modulated  return. 

The  second-order  (power)  spectrum  the  and  fourth-order  spectrum  (trispectrum)  are  special  cases  of 
the  nth-order  spectrum[8,10j.  The  nth-order  spectrum  is  defined  as  the  Fourier  transform  of  the  nth-order 
cumulant.  Therefore,  the  nth-order  spectrum  is  also  called  the  nth-order  cumulant  spectrum.  For  a  sero 
mean  and  stationary  process  the  second-order  and  fourth-order  cumulant  functions  are  as  follows: 


second-order 


fourth-Qrdec 


C7(t)  =  E[x{t)x{t  +  r)|, 


<?«(n»  r3)  =  E\x(t)x(t  +  rt)x(t  +  r2)x(t  +  r3)j 

-J5[x(t)x(t  +  ri))£(x(t  +  r2)x(t  +  r3)) 
-E\x(t)x{t  +  r2)|£|x(t  +  ri)x(t  +  r3)] 
-E\x(t)x(t  +  r3)]£[x(t  +  rj  )x(t  +  r2)|. 
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Since  the  Gaussian  expansion  is  subtracted  from  the  fourth-order  moment,  the  fourth-order  cumulant  is  sero 
for  Gaussian  processes.  Therefore,  as  uoted  by  Rosenblatt  [5j,  additive  Gaussian  noise  will  not  effect  the 
result  of  the  trispectrum. 

This  report  is  concerned  with  a  special  case  of  the  fourth-order  cumulant.  Namely,  when,  =  0,  and 
r2  =  r3.  Therefore,  the  fourth-order  cumulant  reduces  to  the  following: 


C43(r)  =  E\x(t)2x(t  +  r)2)  -  E\x(t)2)E\x(t  +  r)2)  -  2 E[x(t)x(t  +  r)|2 

»  fl3(r)  -  Var[*(0]2  -  2[*(r)|2.  (1.1) 

If  in  addition  r  =  0,  then  C43  is  related  to  the  kurtosis  of  x.  Properties  of  kurtosis  were  discussed  in 
previous  publications|ll,12j. 

The  fourth-order  spectrum  is  therefore, 


/+00 

1 

-OO 


C43(r)e~"Td  r 


R7(r)riutdT 


2x(Var(x(t)|)26(w) 


R(r)2e->urdr. 


(1.2) 


Another  special  case  of  the  fourth-order  cumulant  arises  when  n  »  r3  =  0  or  when  n  —  r3  =  r3  =  r. 
For  this  case, 


C43(r)  =  E[x(t)*x(t  +  r)|  -  3E[x(t)2\E\x(t)x(t  +  r)) 


=  *3(r)  -  3Var(*(t)|R(r). 


When  the  process  is  Gaussian,  C 43(f)  =  0.  Therefore,  its  spectrum  is  also  sero  for  all  frequencies. 

2.  Mixture  Processes 

The  need  for  multidimensional  non-Gaussian  processes  arise  from  the  definition  of  the  nth-order  cumu¬ 
lant  spectra.  Here  non-Gaussian  processes  will  be  modeled  as  non-Gaussian  mixture  processes.  As  alluded 
to  in  the  introduction,  these  mixture  processes  are  more  than  mathematical  conveniences,  because  they  seem 
to  represent  physical  processes. 

A  bivariate  mixture  probability  density  is  defined  as  follows: 

/(*. y)  =  (1  -  A)J/i(z,y)  +  A 2f2(x,y)  +  A(1  -  A)(/,(z)/2(y|z)  +  /2(x)/,(y|z)|,  (2.1) 

2 
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where,  0  <  A  <  1  and  /( x,  y)dxdy  =  1. 

The  bivariate  densities  /i(x,  y)  and  /i(x,  y)  can  be  in  general  any  bivariate  densities.  However,  Gaussian 
densities  are  utilised  in  this  report  in  order  to  obtain  tractable  results.  But,  even  when  /i(x,  y)  and  /2(x,y) 
are  Gaussian  the  mixture  density  /(x,  y)  is  in  general  non-Gaussian.  This  was  demonstrated  for  univariate 
densities  by  Bendat  and  Piersol  [13).  If  /i(x,y)  ^  /2(x,y)  and  A  ^  0,1,  then  /(x,y)  is  non-gaussian.  The 
Gaussian  densities  / i(x,  y)  and  /2(x,  y)  are  defined  below. 

Suppose  the  marginals  are  the  same  in  (2.1),  i.e.,  /i(x)  =  /s(x).  Then, 

/(x,y)  =  (1  -  A)/i  (x,  y)  +  A/2(x,  y).  (2.2) 


As  an  application  of  (2.2),  let  the  mixture  process  be  passed  through  a  hard  clipper.  It  is  desired  to  find 
the  autocorrelation  function  at  the  output.  The  hard  clipper  is  sometimes  used  in  sonar  arrays  and  other 
applications  include  neural  networks[14).  For  the  Gaussian  case  the  result  is  well  known[15,16|.  However, 
for  the  non-Gaussian  case  no  general  result  is  known.  But,  a  general  result  for  the  non-Gaussian  mixture 
density  (2.2)  reduces  to  the  following, 

f?r(r)  =  (1  -  A)-sin_1|pi(r))  +  A-sin~l[p2(r)]. 

If  If 

If,  A  =  0, 1,  then  ft,.(r )  reduces  to  the  Gaussian  case.  Also,  the  autocorrelation  function  at  the  output  of 
more  general  nonlinearities  [17]  can  be  similarly  derived. 

The  form  of  (2.1)  is  more  general  than  the  bivariate  mixture  densities  discussed  in  reference  4.  This 
generality  represents  an  improvement,  since,  (2.1)  does  not  suffer  from  the  disadvantages  of  the  mixture 
densities  of  reference  4.  For  example,  if  x  and  y  are  independent  than  (2.1)  and  (2.2)  reduce  to  a  product 
of  univariate  mixture  densities.  Also,  (2.1)  and  all  special  cases,  e.g.,  (2.2),  always  satisfy  the  strong  mixing 
condition. 

The  form  of  (2.1)  can  also  be  extended  to  a  multivariate  non-Gaussian  mixture  density,  but  this 
generality  is  not  needed  here. 

Let, 


and 


/i(*.  y)  = 
h(x,  y)  = 


l 


2»<Mi  —  P?)l/a 
l  l 


2*<73(i  -  p\yi 2 

be  the  mixture  densities  of  (2.1).  Then  the  autocorrelation  function, 


R{t)  =  f  xyf(x,y)dxdy, 
J  -oo 
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reduces  to  the  following, 


A(r)  =  1(1  -  A)Pi(r)  +  -W»-)!|(l  -  A)crf  +  Acr|)). 


The  spectrum  of  (2.3)  is 


S(u)  =  R(T)e-JUTdr 

J -CO 

f oo  r oo 

=  ((1  -  A)cr2  +  Acr|]  (1  -  A)  /  pi(r)e_-’"rdr  +  A  /  p2(r)e_-’wrdr 

J  —  oo  J  —  oo 


On  the  other  hand,  the  fourth-order  moment  of  (2.1), 


reduces  to  the  following  result, 


/oo 

x2y2f(x,y)dxdy, 

-  OO 

ftz(r)  =  ((1  -  a)<ti  +  ^1? 

+2[(1  -  A)p?(r)  +  Apf(r))((l  -  A)<rf  +  Acr^J 
+A(l-A)|p?(r)<rI-p|(r)rr?|(<rI-<r?). 


Since  R2(r)  is  clearly  different  from  a  linear  combination  of  R(r)2  a  non-Gaussian  mixture  process  can 
be  differentiated  from  a  Gaussian  process. 

From  the  definition  of  the  fourth-order  cumulant  given  in  (1.1), 

C«(r)  =  R2(r)  -  ((1  -  A)<r2  +  A<r||2  -  2R(r)2.  (2.6) 

As  expected,  if  A  =  1,0,  (the  Gaussian  case),  then  C42(r)  =  0.  However,  for  the  non-Gaussian  case, 
when  <Ti  =  cr2  =  1, 

C43(r)  =  2A(1  -  A)(pt(r)  -  p2(r)|2.  (2.7) 

The  fourth-order  spectrum  of  (2.7)  is, 


C42(w)  =  2A(1  -  A)  f  |pi(r)  -  p2(r))2e  ,UTdr. 

J  -OO 


EXAMPLE  1. 


Let,  Pi(r)  =  e  and  p2(r)  =  e  in  equation  (2.8),  then  for  A  /  0, 1, 


,  _  ^43 (w)  _  4qi  4a2 _ 4(ai  +  a2) 

42'W'  ~  2A(1  -  A)  ~  w2  +  (2ai)2  w2  +  (2q2)2  u>2  +  (a,  +  a2)2 ' 


Equation  (2.9)  is  plotted  for  positive  frequencies  in  figure  1  for  three  cases:  (c«i  =  10,  =  l),(ai  = 

10,  =  2),  and  (an  —  10,  aj  =  9).  The  last  case  is  close  to  sero  for  all  frequencies  indicating  that  the 

process  is  almost  Gaussian.  The  other  two  cases  show  that  the  fourth-order  spectrum  can  differentiate 
between  Gaussian  and  non-Gaussian  processes.  It  is  also  apparent  from  the  figure  that  the  fourth-order 
spectrum  is  negative  for  some  frequencies. 

EXAMPLE  2. 

Let,  pi(r)  =  e-'Mlf|cos(a»0r),  and,  p3(r)  =  e_a,lrlcos(u;or),  in  (2.8),  then,  for  A  ^  0, 1, 

n'  (a,!  =  g*2H  _  2ai  .  _ _ 

2A(1  —  A)  u 2  +  (2ai)3  (w  —  2uio)^  +  (2c«i)2  (ui  +  2wo)2  +  (2ai)3 

2a2  aj  aj 

w2  +  (2 aj)2  (ui  —  2wo)2  +  (2aa)2  (w  4-  2u?o)2  +  (2a3)2 

2(ai  -f  a3) _ ai  +  a2 _ at  +  a2 _  . 

a;2  +  (ai  +  a2)2  (w  —  2wo)2  +  (ai  +  a3)2  (t*>  +  2u>o)2  +  (ai  +  c«2)2 

Equation  (2.10)  is  plotted  in  figure  2  for  the  same  conditions  as  given  in  example  1. 

The  fourth-order  moment,  for  the  other  special  case  mentioned  in  the  introduction,  is 
C43(r)  =  3A(1  -  A)[p3(r)<r2(cx2  -  <r|)  +  pi(r)<r£(o2  -  a2)]. 

If  tri  =  <T2.  then  C43(r)  =  0.  This  case  will  not  be  considered  further. 

3.  Sinusoidal  Mixture  Processes 

A  problem  that  occurs  in  underwater  acoustics  can  be  stated  as  follows:  given  s(t),  -T<t<T,  determine 
if  s(t)  represents  a  sum  of  two  sinusoids  or  a  mixture  of  two  sinusoids.  A  more  complicated  problem  arises 
from  creeping  wave  backscatter.  This  case  is  usually  discussed  in  the  context  of  target  classification  from 
resonances)  18].  Here,  however,  the  backscattered  waveform  is  idealised  as  a  mixture  of  two  sinusoids  in 
order  to  show  the  classification  potential  of  fourth-order  spectra.  Whereas,  both  a  sum  and  a  mixture  of  two 
sinusoids  will  have  fourth-order  spectra,  only  the  mixture  represents  the  idealised  backscattered  waveform. 
Therefore,  the  fourth-order  spectrum  gives  additional  information  about  the  physical  phenomenon  that  can 
be  utilised  in  target  classification. 

Nevertheless,  it  is  interesting  that  the  idealised  problem  can  be  solved  by  fourth-  order  spectra.  On  the 
other  hand,  second-order  spectra  cannot  differentiate  between  a  sum  and  a  mixture,  since  both  could  have 
the  same  spectra.  However,  there  may  be  amplitude  differences,  but  this  information  is  usually  not  enough 
to  make  a  determination. 
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It  will  be  assumed  that  the  data  are  ergodic  so  that  the  time  autocorrelation  function  !R(r)  and  the 
time  fourth-order  function  8C2 (r)  of  the  random  process  equals  the  corresponding  statistical  autocorrelation 
function  and  statistical  fourth-order  function,  respectively)  19|. 

a)  Sum  of  Sinusoids 

Let, 

a(t)  =  ^sin(cj1t  +  <£i)  +  ^=sin(w2t  +  ^2),  (3.1) 

where,  wj  5^  a i2  ,  —  T  <  t  <  T,  and  the  parameters  di  and  fa  are  random  variables  uniformly  distributed 
between  0  and  2ir. 


The  time  autocorrelation  function, 


reduces  to, 


*w  * r'it  2<rrwi S-rl +  **• 

2  2 

R(r)  =  T'co«(wir)  +  — cos(w2r)  =  R(r). 

4  4 


Its  corresponding  spectrum  is  given  in  (3.3), 

2ko[  £(w  -  wi)  4-  S(u>  +  o>i)  2ira%  5(w  -  w2)  +  +  w2) 

~  4  2  +  4  2 


Similarly,  the  time  fourth-order  function, 


*lW  ‘  />»  lirTwl  I-r'l  + 


for  (3.1),  reduces  to, 


»2(r)  =  a~+  2°1^>2  +  +  ^cos(2wir)  4-  j~cos(2w2r)  +  ^±pcos(wir)cos(u>2r)  =  R2(t). 


Notice  that  the  spectrum  of  R2(r)  will  have  5  positive  frequencies.  However,  the  fourth-order  cumulant 
is  desired.  It  reduces  to, 


C*i(r)  =  ^  "  j^cos^r)  -  ^cos(2w2r), 

and  its  corresponding  spectrum  follows  as, 

C  (qi)  -  ^’r(ai  +  aj)  j(^|  2^^  S(w  -  2wi)  +  6(w +2u>i)  2ira\  6(a/  -  2w2)  +  S(w  +  2a>2) 
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Therefore,  the  fourth-order  cumulant  spectrum  for  (3.1)  contains  only  3  positive  frequencies  located 
at  0,2wi,  and  2w2.  The  frequency  locations  for  the  fourth-order  cumulant  of  the  mixture  process  will  be 
different.  This  represents  a  means  to  differentiate  a  sum  from  a  mixture. 

b)  Mixture  of  Sinusoids 


Let, 

s(t)  =  aisin(wtt  +  01)hl(t)  +  a2sin(tu2t  +  <£2)h2(t),  (3.6) 

represent  the  sinusoidal  mixture  process  over  the  interval,  —Tl  <  t  <  Tv,  where  Tl  >  0,  and  Tv  >  0.  The 
functions  hi(t)  and  h2(t)  are  defined  in  the  following  way, 


hl(t)  =  u(t  +  TL)-u(t-Tl), 

and 

h2(t)  =  u(t  —  Ti)  -  u(t  -  Tv), 

where  u()  is  the  unit  step  function.  These  functions  also  have  the  property  that  h?()  =  /»,(),  i=  1,2,  which 
will  be  used  later. 


The  time  autocorrelation  function  of  (3.6), 

*W  *  r  Jr-  TWli-ZW  C'l  +  '>"■ 


reduces  to,  for  r  >  0, 


A  similar  expression  is  obtained  for  r  <  0.  These  expressions  converge  to  the  following, 


2  2 

R(r)  =  (1  -  A)ycos(w2r)  +  A^-cos(wir), 

where,  — *  A.  If  A  =  A,  then  (3.7)  and  (3.2)  will  have  identical  spectra. 

The  fourth-order  moment  can  be  shown  to  converge  to  the  following  result, 

_  ,  .  o}A  +  a2(l-A)  ai(l-A)  ,  ,  of  A 

»2(r)  =  - - - -  +  - -cos(2w2r)  +  -A-cos(2«1r). 


(3.7) 


(3.8) 
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Notice  that  only  3  positive  frequencies  are  present  in  (3.8).  The  fourth-order  cumulant  is  given  by, 


a} A2  aj(l  -  A)2 


C*2(r)  =  JA(1  -  A)(a2  -  a2)2  -  -  4 

-cos(2wir)  - 


S  <»j(l  ~  A)(l  —  2A)  _  ln  ,v  ,  2 


cos(2w2r)  -  a2a|A(l  -  A)cos(wir)co8(a)2r). 


If,  ai  =  a2  =  1,  and  A  =  A,  th«  fourth-order  cumulant  reduces  to, 


<?42(r)  =  -  jcos(w1r)co8(w2r), 


(3.9) 


and  its  corresponding  spectrum, 


^  ,  .  2ir  2*  {(<■/  — (<«>2  —  <i>i)j  +  6(u  +  (fa/?  ~  u>i))  2x  5[w  —  (w2  4-  wj)|  +  6[u>  +  (ui2  +  wi)l 

C*7'w>  -  -—b(ui)  -  -  2  1  2  ’ 

(3.10) 

is  clearly  different  from  the  spectrum  in  (3.5),  for  the  same  conditions. 


EXAMPLE  3. 


(  Sum  of  Sinusoids) 

Figure  3  represents  the  spectrum  (in  db)  of  (3.2)  for  at  =  a2  =  1,  and  A  =  A.  The  corresponding 
fourth-order  cumulant  spectrum  (for  a  linear  scale)  is  shown  in  figure  4. 

(  Mixture  of  Sinusoids) 

The  spectrum  of  (3.7)  is  shown  in  figure  5,  for  «i  =  a2  =  1,  and  A  =  A.  In  this  case  the  spectra  of 
figures  3  and  5  are  identical.  However,  the  fourth-order  cumulant  spectrum  (3.10)  as  shown  in  figure  6  is 
clearly  different  from  figure  4.  Therefore,  the  sinusoidal  mixture  process  can  be  differentiated  from  a  sum  of 
sinusoids,  based  on  the  fourth-order  cumulant  spectrum. 

4.  Modulated  Processes 


Modulated  processes  arise  in  communications  as  well  as  sonar  and  radar  applications.  Here  the  following 
general  modulated  problem  is  treated.  Let, 


x(t)  =  a(t)z(t),  (4.1) 

be  the  modulated  process.  Where  a(t)  and  s(t)  are  both  sero  mean  and  mutually  independent  processes. 
Equation  (4.1)  is  in  general  a  non-Gaussian  process  even  if  both  a(t)  and  s(t)  are  Gaussian  processes.  Its 
autocorrelation  function, 

Rx(t)  =  E[x(t)x(t  +  r)|  =  E\a(t)a(t  +  r)\E[z(t)z(t  +  r)|  =  f?„(r)ff,(r),  (4.2) 
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is  a  product  of  two  autocorrelation  functions.  Therefore,  its  spectrum  will  be  a  convolution  of  the  two 
spectra.  If  s(t)  is  the  information  bearing  component  then  the  modulating  component  a(t)  tends  to  interfere 
with  the  reception  of  s(t). 

Let  a(t)  be  a  Gaussian  process.  Then  the  fourth-order  cumulant  of  equation  (4.1)  reduces  to 

C«2(r)  =  Var2(a(t)|  [E[z(t)3z(t  +  r)2|  -  Var2[s(t)|j  +  2Ra(r)2  [£(z(e)2z(t  +  r)2|  -  fl,(r)2] .  (4.3) 

Notice  that  (4.3)  is  not  a  simple  product  of  autocorrelation  functions.  Therefore,  its  spectrum  will  be 
different  from  the  spectrum  of  (4.2). 

Now  let  z(t)  be  Gaussian  also.  Then  (4.3)  reduces  to, 

Cu(r)  =  2Var2[a(t))J?1(r)2  +  2Var2[*(t)|,R„(r)2  +  R<,(r)2  R,(r)2 ,  (4.4) 

which  can  be  seen  to  separate  the  two  spectra  into  a  sum  and  a  product  rather  than  just  a  product. 

EXAMPLE  4. 

Let, 

x(t)  =  a(t)cos(w0t  +  4), 

where,  ^  is  a  random  parameter  uniformly  distributed  between  0  and  2x.  Then, 


If, 


R.,{r)  =  e-a|r|, 


then  the  spectrum  of  Rx(r)  is, 


_  .  .  a  a 

5,(w)  =  t ;  + 


(w  -  w0)2  +  a2  (w  +  wq)2  +  a2 


(4.5) 


If  a  -♦  oo  ,  then  Sx(ui)  — *  0.  This  means  that  the  sinusoidal  frequency  wo  will  not  be  observed  in  the 
spectrum  of  Rx(r). 


However,  the  fourth-order  cumulant, 

C«(r)  =  ^fl,.(r)2  +  — — -~-^cos(2w0r), 

has  the  following  spectrum, 


Cu(u>)  = 


a 

w2  -I-  (2a)2 


2s-Var2[o(t)j  5(w  -  2w„)  +  £(w  +  2wo) 
+  8  2 


(4.6) 


Now,  the  sinusoidal  frequency  2w„  is  seen  in  (4.6)  even  as  a  — »  oo. 
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(Simulation:  Amplitude  Modulated  Sinusoid) 

The  spectrum  (in  db)  of  a  sinusoid,  amplitude  modulated  by  white  noise,  is  shown  in  figure  7.  The 
sinusoidal  frequency  is  not  discernable  in  this  figure  because  of  the  interference  of  the  modulation  as  predicted 
in  (4.5).  However,  the  fourth-order  spectrum  in  figure  8  (plotted  on  a  linear  scale)  reveals  the  underlying 
sinusoid  but  at  twice  its  frequency.  This  result  is  predicted  in  (4.6). 

EXAMPLE  5.  (  SONAR/RADAR  APPLICATION  ) 

Let, 

y(t)  =  £  M t-nTr)e*“'‘+#«', 

n=  —  k 

-oo  <  t  <  oo,  be  the  pulse  coded  transmitted  signal  of  length,  2kTp  +  T.  The  pulses  are  defined  as  follows, 

h(t  -  nTr)  =  u(t  -  nT„  +  j)  -  u(t  -  n7),  -  y), 

where  u()  is  the  unit  step  function.  The  parameters,  T  and  Tp,  for  T  <  T„,  represent  the  pulse  width  and 
pulse  repetition  interval,  respectively. 

The  received  signal  is  assumed  to  be  of  the  form, 

x'(t)  =  o(t  -  ~)  Y'  h(t  -  nTP  -  r«)eJl,w*+w'*<"*-r»,+9'‘+*l, 

n=-fc 

where,  Tp  is  the  range  of  the  target,  is  the  Doppler  shift  radian  frequency  associated  with  each 
transmitted  radian  frequency  u„,  8n  is  the  phase  of  the  nth  transmitted  pulse,  and  <f>  is  a  random  phase  angle 
uniformly  distributed  between  0  and  2tt.  The  real  stochastic  Gaussian  modulating  function  a(t)  represents 
a  time  fluctuating  target(20j.  This  model  usually  applies  for  transmitted  signals  of  long  length. 

Let, 

y'(f)  *  £  h(t-nTP-rR)eJ^t~T^+t'K 

f»=  —  fc 

be  the  adjusted  transmitted  signal.  Where  Vp  is  a  parameter  that  is  adjusted  in  order  to  search  for  the  true 
range.  For  simplicity,  it  will  be  assumed  that  T'r  =  Tr  ,  i.e.,  the  range  is  known. 

The  envelope  of  the  received  signal  is  defined  as  follows, 

z(t)  =  z'(t)y'*(t), 


where  the  asterisk  is  the  complex  conjugate. 
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Therefore,  since  the  pulses  are  disjoint,  i.e., 

the  envelope  of  the  received  signal  is, 


i(t)  =  a(t  - 


h{t  -  nTv  -  rR)e^w''‘,‘-r",+^. 


(4.7) 


The  objective  of  this  example  is  to  compare  the  autocorrelation  function  of  (4.7)  and  its  corresponding 
spectrum  with  the  fourth-order  cumulant  function  of  (4.7)  and  its  corresponding  spectrum. 

The  autocorrelation  function  is  given  by  the  following  expression, 

k  k 

£|i(t1)x*(t2))  =  R,.(ta-t1)  Y,  — .«*»>,  (4.8) 

ill  =  -k»i  = -It 


where, 


R«(h  -  tt)  =  fi[a(*i  -  ^)o(t,  -  ^)), 


is  a  stationary  autocorrelation  function. 


Therefore,  since  (4.7)  is  not  in  general  a  stationary  process,  the  spectrum  of  (4.8)  is  defined  as  a  two- 
dimensional  Fourier  transform  of  the  autocorrelation  function, 

S(TRtoi,u,)=  r  r  E\x(ti)x'[t?)]e~liuitl+W,t,idtldt2.  (4.9) 

J-OO  J -oo 


To  simplify  the  evaluation  of  (4.9)  the  modulating  function  will  be  assumed  to  be  white  noise,  i.e.,  /?„(t2  - 
)  =  6(t2  -  1 1).  Therefore,  the  spectrum  of  (4.8)  reduces  to 


S(Tb,Wi,w2)=T«-j(,-‘+“j't" 


sin[(wi  +wa)£|sin[(m  +  u>2)(fc  +  ^)r,,[ 
[(wi+u>2)£|  sin[(w!  +w2)^] 


(4.10) 


The  spectrum  (4.10)  is  concentrated  along  the  line,  wi  =  -w2  with  the  same  value  ,  T(2k+1),  for  all 
frequencies.  Therefore,  if  =  —  w2,  then,  S(TR,w)  =  T(2k  4-  1),  Vu>.  This  clearly  shows  the  interfering 
effect  of  the  modulating  white  noise. 

The  fourth-order  cumulant  function  reduces  to  the  following  result, 

k  k 

<?42(t„t2)  =  Var2(a|  £  £  A(‘i  *  ntTp  -  TR)h{ti  -  n7Tp  -  (4.11) 

HlS-fc  f»3  =  —  fc 
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where  the  Gaussian  assumption  for  the  modulating  function  was  used  to  obtain  (4.11). 

Similarly,  the  fourth-order  cumulant  spectrum  is  defined  as  a  two-dimensional  Fourier  transform, 

/OO  A  oo 

/  (4.12) 

■oo  J  -oo 

which  reduces  to, 


C*7(Tn,wuui)  =  Var2|a]T2  ^ 

ft  1  =  —k  fljS5“fc 

sin [(uii  -  2a;ni,t)^l  sin[(ai2  4-  2aJn,.t)|] 


[  (ail  -  2a i„„j)f  (w2  +  wn,rf)f  J  ^  ^ 

Again,  the  spectrum  is  concentrated  along  the  line,  o/i  =  —013,  but  now  with  peaks  at  twice  the  Doppler 
shift  frequencies.  To  see  this  more  clearly,  let  oil  =  -oi2  =  ui  ,  and  uin,i  =  w,i,Vn,  i.e.,  all  the  Doppler  shift 
frequencies  are  the  same,  then 


C«(r*,w)  =  Var2[a)r2(2Jfc+l)2 


sin((a;  -  2a>rf)|| 


(w  -  2wj)^ 


(4.14) 


It  is  clear  that  (4.13)  and  (4.14)  are  not  interfered  with  by  the  modulating  fun- ‘ion,  except  for  the 
constant  representing  its  variance.  Therefore,  the  fourth-order  cumulant  spectrum  of  the  received  signal 
will  be  discernable.  Whereas,  the  spectrum  of  (4.8)  may  not  be  discemable  depending  on  the  modulating 
function  a(t)  as  shown  in  example  4. 

5.  Conclusions 

It  has  been  shown  that  for  non-Gaussian  processes,  modeled  as  mixture  and  modulated  processes, 
higher-order  spectra  can  extract  information  that  second-order  spectra  cannot.  Namely,  the  fourth-order 
cumulant  spectrum  can  differentiate  between  purely  Gaussian  processes  and  mixture  non-Gaussian  processes. 
This  was  demonstrated  by  several  examples.  Moreover,  the  fourth-order  cumulant  spectra  can  differentiate 
between  a  sum  of  sinusoids  and  a  mixture  of  sinusoids.  This  was  demonstrated  by  simulation. 

An  important  practical  application  of  fourth-order  spectra  showed  that  target  information,  from  an 
active  sonar  return,  could  be  extracted  .  Specifically,  the  fourth-order  spectrum  extracted  Doppler  shift 
frequencies  from  an  amplitude  modulated  coded  pulse  train  return,  whereas,  the  spectrum  was  unable  to  do 
so. 
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Spectrum:  Oj  =  10. 
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FIGURE  7. 
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